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Analytical  and  experimental  investigations  of  the  ignition  and 
combustion  of  hydrogen-oxygen-argon  mixtures  are  presented.  A  one¬ 
dimensional  kinetics  program  with  generalized  chemistry  and  provisions 
for  mass  addition,  momentum  addition,  heat  loss,  mixing,  shock  waves, 
and  a  rate  screening  option  has  been  developed  and  used  to  analyze  the 
effect  of  free  radical  additives  on  the  ignition  delay  time  in  hydrogen-oxygen 
mixtures.  This  computer  program  has  also  been  used  to  reduce  shock  tube 
measurements  of  hydrogen-oxygen  ignition  delay  and  OH  S  -  n  emission. 
The  experimental  results  indicate  that  the  reaction  O  +  H  (+  M)  -*  OH*  (+  M) 
is  responsible  for  producing  OH  (^~b)  during  the  induction  period;  however, 
this  mechanism,  when  input  into  the  computer  progiam,  was  not  sufficient 
to  qualitatively  reproduce  OH  emission  intensity  profiles  obtained  experi¬ 
mentally.  Work  aimed  at  determining  the  mechanism  for  OH  emission  through 
out  the  reaction  zone  is  continuing . 
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Combustion  in  a  supersonic  stream  has  become  one  of  the  most  important 
new  approaches  to  achieving  hypersonic  flight  with  air  breathing  engine  systems. 
In  the  past  few  years,  a  considerable  amount  of  effort  has  been  expended  in 
theoretical  analysis  of  the  combustion  process  and  in  the  design  of  hardware. 

A  number  of  fundamental  questions  still  remain  unresolved,  however.  These 
questions,  pertaining  to  the  mechanism  of  ignition  and  reaction  in  a  rapidly 
mixing  fuel-air  system,  the  establishment  of  suitable  reaction  rate  expressions 
under  conditions  where  reaction  and  mixing  occur  at  competitive  rates ,  and  the 
effect  of  additives  on  the  chemically  reacting  system,  are  still  essential  and 
yet  unresolved  problems. 

This  report  describes  the  work  conducted  under  AFOSR  Contract  No.  F44620- 
68-C-0069  to  gain  insight  into  these  fundamental  processes.  This  work 
is  being  continued  under  AFOSR  Contract  No.  F44620-70-C-0061 . 

Analytical  Investigations  of  ignition  and  combustion  and  the  effect  of 
additives  on  ignition  delay  were  carried  out  for  the  hydrogen-oxygen  system 
using  an  existing  one-dimensional  chemical  kinetics  computer  program 
significantly  modified  to  include  generalized  chemistry,  mass  addition, 
momentum  addition,  heat  loss,  mixing,  shock  waves,  and  the  capability  of 
conducting  a  chemical  rate  screening  analysis. 

An  experimental  program  was  conducted  in  order  to  better  understand 
the  details  of  the  complex  chemical  processes  occurring  during  ignition 
and  combustion.  The  portion  of  the  experimental  program  completed  to  date 
consists  of  measurement  of  OH  -  n  emission  for  both  the  ignition  zone 
and  the  entire  combustion  zone.  These  measurements  were  used  to  determine 
the  ignition  delay  times  and  are  currently  being  examined  in  conjunction  with 
the  computer  program  to  deduce  the  mechanism  that  produces  excited  OH. 
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COMPUTER  PROGRAM  DEVELOPMENT 

At  the  outset  of  this  program,  it  was  decided  that  the  development  of 
a  complex  computer  program  which  treated  the  fluid  mechanical  details  of 
turbulent  mixing  and  diffusion  in  as  exact  a  manner  as  is  possible  by  the 
current  state-of-the-art  would  be  both  duplicative  (several  programs  of  this 
type  are  currently  in  existence,  for  example,  Reference  1),  and 
beyond  the  scope  of  our  particular  line  of  investigation;  therefore,  it  was 
decided  to  develop  an  analysis  which  treated  the  chemical  kinetics  exactly 
and  had  provisions  for  mass  transfer,  heat  tranfer  and  momentum  addition  in  a 
simplified  manner.  Above  all,  this  program  would  be  simple  to  use,  require  a 
minimum  of  computer  time,  and  be  easily  adaptable  to  handle  diverse  problems. 
It  was  envisioned  that  this  tool  would  find  its  maximum  usefulness  during  the 
development  and  testing  stage  of  actual  scramjet  hardware.  These  objectives 
have  been  achieved. 

The  computer  program  utilizes  the  JANNAF  (ICRPG)  One-Dimensional 
Kinetics  Program  as  a  framework  to  which  significant  modifications  were 
added,  greatly  enhancing  its  capabilities  as  a  diagnostic  tool.  The  major 
features  of  the  program,  as  it  presently  exists,  are  briefly  enumerated  below, 
and  a  detailed  description  of  the  program  mechanics  can  be  found  in  Appendix  A. 

1 .  The  method  of  solution  has  been  modified  to  allow  the  integration 
of  the  conservation  equations  for  the  chemical  system  in  a  form 
such  that  the  chemistry  is  generalized.  This  generalized  chemistry 
feature  allows  the  user  to  easily  input  up  to  40  species  and  158 
reactions  in  symbolic  form,  including  specified  third  body 
efficiencies  for  each  reaction,  if  necessary.  The  JANNAF  thermo¬ 
chemical  data  is  available  as  part  of  the  program  in  a  form  which 
allows  for  future  expansion  directly  from  the  JANNAF  format. 
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2.  The  fully  Implicit  numerical  integration  method  used  in  the  program 
makes  it  the  most  efficient  program  of  its  type,  particularly  when 
handling  chemical  systems  near  equilibrium  or  in  an  ignition 
region.  It  has  been  shown  that  explicit  methods  of  numerical 
Integration  are  unstable  when  applied  to  relaxation  equations 
unless  the  integration  step  size  is  of  the  order  of  the  character¬ 
istic  relaxation  distance.  Since  in  the  near  equilibrium  flow 
regime  the  characteristic  relaxation  distance  is  typically  many 
orders  of  magnitude  smaller  than  characteristic  physical 
dimensions  of  the  system  of  interest,  the  use  of  explicit  methods 
to  integrate  relaxation  equations  often  results  in  excessively  long 
computation  times.  An  Implicit  integration  method  which  is 
inherently  stable  in  all  flow  situations  (whether  near  equilibrium 
or  frozen)  is  therefore  used  by  the  computer  program.  With  this 
method,  step  sizes  which  are  of  the  order  of  the  physical 
dimensions  of  the  system  of  interest  can  be  used,  reducing  the 
computation  time  and  cost  per  case  several  orders  of  magnitude 
when  compared  with  conventional  explicit  integration  methods . 

3.  Mass,  momentum,  and  energy  addition  rates  can  be  specified  as 
part  of  the  initial  Input  conditions . 

4 .  Finite  mixing  times  are  accounted  for  by  holding  back  a  portion  of 
the  initial  fuel  from  reaction,  until  after  a  finite  "mixing  time"  or, 
when  a  specified  axial  position  is  reached. 

5.  If  shock  waves  exist  at  known  locations  in  the  flow,  the  program 
will  consider  them  exactly  if  the  program  input  contains  either 
the  downstream  pressure,  velocity  or  the  shock  angle. 

6.  The  program  has  also  been  modified  so  as  to  make  available  a 
unique  analysis  capability  for  complex  reacting  systems.  At 
the  desire  of  the  user  the  program  will  print  out,  in  increasing 
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order,  the  quantitative  contribution  of  each  reaction  to  the 


production  of  any  of  the  species  taking  part  in  the  reacting 
system.  The  program  will  then  eliminate  from  the  reaction 
set  those  reactions  whose  relative  contribution  to  the 
production  of  any  specie  falls  below  a  preset  level.  The 
integration  step  is  then  recomputed  with  the  reduced 
reaction  set  and  all  parameters  printed  out  along  with  the 
original  computation  for  comparison. 


It  should  be  emphasized  that  one  of  the  most  significant  features  of 
this  program  is  its  operating  efficiency.  While  a  number  of  options  are 
available  which  makes  the  program  an  extremely  flexible  tool,  the  required 
information  needed  to  automate  any  particular  option  is  straightforward  and 
can  be  input  easily.  In  addition,  the  program's  very  low  operating  costs 
per  computer  run  make  it  ideal  for  use  when  a  large  number  of  cases  must 
be  run,  for  example,  during  the  development  and  testing  stage  of  actual 
hardware . 


Appmxs  EffKXWEWfss-STjpiaa 

A  theoretical  study  of  the  effect  of  additives  on  the  ignition  delay  of 
the  hydrogen-oxygen  reaction  has  been  presented  in  Reference  2 . 

Since  that  work  was  presented,  the  program  was  further  modified  by  the 
addition  of  a  reaction  screening  option  which  allowed  an  even  more  detailed 
investigation  of  the  entire  reaction  process.  Using  this  option,  it  is  possible 
to  print  out  the  specific  contribution  of  each  reaction  in  the  chemical  system 
toward  the  production  of  each  of  the  species  that  exist  in  the  system.  This 
data  immediately  pinpoints  the  dominant  reaction  or  reactions  at  any  time 
during  the  chemical  process .  This  information  must  be  available  if  reaction 
rate  constants,  or  reaction  mechanisms  are  to  be  obtained  from  experimental 
measurements  made  during  a  chemical  process  involving  a  multitude  of 
individual  reactions . 
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Figure  1  which  plots  the  rate  of  production  or  destruction  of  hydrogen 
atoms  for  each  reaction  versus  time  for  given  initial  conditions,  illustrates 
the  ease  with  which  the  details  of  the  chemical  process  can  be  seen  by  this 
type  of  presentation.  From  this  figure  it  can  be  seen  that  the  production  of 
the  dominant  chain  carrier  (H  atoms)  during  the  induction  period  of  the 
hydrogen-oxygen  reaction ,  is  controlled  by  the  three  reactions 

H2  +  OH  =  H20  +  H 
H2  +  O  =  OH  +  H 
02  +  H  =  OH  +  O 

In  addition,  the  transition  of  the  hydrogen  dissociation  reaction 
H„  -♦  2H,  to  a  recombination  reaction  H  +  H  -♦  H„ ,  can  be  seen  to  occur 

u  4 

at  approximately  65  microseconds  with  a  similar  reversal  occurring  for  the 
reaction  C>2  +  H  -♦  HOz  at  125  microseconds.  The  chemical  mechanism  which 
was  used  for  this  calculation  is  given  in  Table  1. 

It  is  in  the  region  past  "ignition'*  when  the  contribution  of  several 
reactions  to  the  production  of  H  atoms  are  the  same  order  of  magnitude,  that 
the  screening  option  of  the  program  is  particularly  useful.  If  another  run  is 
made  with  the  same  initial  conditions  but  with  a  modified  reaction  rate 
parameter  for  a  particular  reaction,  the  effect  of  this  perturbed  rate  on  the 
total  reaction  picture  will  become  readily  apparent.  This  option,  therefore, 
can  be  used  to  set  limits  within  which  rate  parameters  must  be  known  in  order 
to  have  a  particular  degree  of  confidence  in  the  computed  results .  The  program 
can  also  be  used  to  establish  physical  conditions  for  which  reactions  whose 
rate  parameters  are  poorly  known  and  possibly  difficult  or  impossible  to  measure 
in  an  isolated  system,  are  dominant  and  control  the  variation  of  a  readily 
measured  physical  property  of  the  complex  system,  e.g. ,  another  specie 
concentration,  T,  p  ,  etc.  Experiments  can  then  be  conducted  at  these 
conditions  to  yield  data  which  can  be  directly  and  quantitatively  linked  to 
the  unknown  reaction  rate. 
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The  screening  option  was  used  to  investigate  the  details  underlying  the 

effectiveness  of  free  radical  producing  additives  in  reducing  ignition  delay. 

Figure  2  shows  the  results  of  2  calculations  for  identical  initial  conditions 
—8 

except  for  6x10  mole  fraction  of  H  atoms  added  for  the  second  calculation. 

As  can  be  seen,  the  effect  of  this  "additive"  on  the  ignition  and  combustion 

process  is  to  shift  the  production  curve  toward  earlier  times.  In  fact,  the 

magnitude  of  this  time  shift  is  just  the  time  that  it  takes  the  reactions 

—8 

alone  to  produce  a  6  x  10  mole  fraction  of  H  atoms.  While  this  result 
may  appear  to  be  obvious ,  it  is  significant  to  observe  that  no  shifts  in 
reaction  mechanism,  or  in  the  dominance  of  any  particular  reaction,  have 
been  effected  by  the  inclusion  of  the  additive.  Therefore,  by  making  one 
analytical  calculation,  it  is  possible  to  know  exactly  what  concentration 
of  additive  must  be  added  initially  in  order  to  achieve  a  particular  reduction 
in  ignition  delay  time. 

The  above  conclusion  will  not,  of  course,  be  valid  for  very  large 
additive  concentrations ,  and  the  limiting  concentration  below  which  these 
results  are  applicable  is  currently  being  determined. 

Experimental  verification  of  this  phenomena  is  planned  for  future 
investigation . 

EXPERIMENTAL  STUDIES 

The  experimental  portion  of  this  program  was  undertaken  to  provide 
data  to  verify  the  predictions  of  the  analytical  model,  and  to  develop 
empirical  relations  which  could  be  used  to  specify  the  mixing  times  or  lengths 
required  as  input  by  the  computer  program.  The  original  plan  for  this  test 
program  involved  monitoring  the  concentrations  of  several  chemical  species 
during  the  ignition  and  combustion  process.  For  ignition  delay  studies,  a 
measurement  of  the  concentration  of  OH  would  be  sufficient  to  determine  the 
characteristics  of  the  combustion  process  for  hydrogen-oxygen.  However, 
for  fuels  other  than  hydrogen,  and  in  particular,  the  planned  mixing  zone 
experiments,  it  was  desirable  to  monitor  simultaneously  the  concentration  of 
two  or  more  species.  Spectroscopic  absorption  measurements  have  been  widely 
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used  to  quantitatively  measure  species  concentration  in  reacting  gas  mixtures. 
However,  because  of  the  difficulty  involved  in  making  simultaneous  absorption 
measurements  at  different  wave  lengths,  and  because  the  experimental  simplicity 
of  emission  measurements  would  enhance  their  use  for  examining  a  multitude  of 
combustion  processes,  it  was  decided  to  use  a  measurement  of  the  emitted  band 
radiation  as  an  indication  of  the  progress  of  the  chemical  reactions. 

A  description  of  the  shock  tube  facility  in  which  the  experiments  were 
conducted  is  given  in  Reference  2  and  is  included  in  this  report  for  completeness. 

In  operation,  test  samples  are  prepared  mano metrically  and  fed  into  the 
driven  section  of  the  tube  through  a  port  in  the  test  section.  During  the  check-out 
phase  of  the  facility,  it  was  found  that  the  burst  pressure  of  the  shock  tube 
diaphragm  could  be  held  to  within  1%  when  a  large  number  of  diaphragms  were 
identically  scribed  on  a  milling  machine.  Throughout  the  experimental  study, 
the  procedure  adopted  was  to  select  a  diaphragm  material,  thickness  and  scribe 
depth  which  would  have  optimum  opening  characteristics  at  a  desired  burst 
pressure.  A  range  of  shock  speeds  were  then  obtained  by  replacing  a  portion 
of  the  hydrogen  driver  gas  with  nitrogen.  Varying  the  molecular  weight  of  the 
driver  gas  in  this  manner,  provided  accurate  control  of  the  incidence  shock 
velocity  throughout  a  1000°K  variation  in  temperature  behind  the  shock  wave. 

The  shock  velocity  was  determined  by  displaying  amplified  outputs  from 
platinum  thin  film  gauges  (located  at  intervals  along  the  tube)  on  a  rastered 
sweep  displayed  on  a  Tectronix  type  556  oscilloscope.  The  raster  was  intensity 
modulated  by  a  Tectronix  type  184  time  mark  generator,  giving  shock  wave 
transit  times  between  gauges  that  were  measurable  to  within  one  microsecond. 

After  a  preset  delay  time  provided  by  the  internal  delay  generator  on  the  556,  a 
Tectronix  type  533A  oscilloscope  was  triggered  which  displayed  the  output  of 
an  EMI  9601B  photomultiplier  tube  (S— 1 1  (O)  extended  range  photocathode 
which  monitored  the  T>  -  II  band  radiation  from  the  OH  molecule  through  a 
Hilger-Engis  Model  600  monochrometer-spectrometer.  Coincident  with  the 
triggering  of  the  type  533A  oscilloscope,  a  signal  is  imposed  on  the  raster 
sweep  thereby  providing  for  timing  correlation  between  the  shock  position  and 
the  ignition  trace.  Sample  oscilloscope  outputs  are  shown  in  Figures  3,4. 
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The  ignition  delay  is  computed  by  measuring  the  time  between  the  last 
heat  tranfer  signal  and  the  initiation  of  the  type  533A  trace.  To  this  time  is 
added  the  time  interval  between  the  start  of  the  type  533A  trace  and  the  first 
appearance  of  an  emission  signal,  measured  directly  from  the  photomultiplier 
output  trace.  The  time  it  takes  the  shock  wave  to  travel  from  the  last  heat 
transfer  gauge  to  the  observation  station  is  calculated  from  the  measured  shock 
wave  velocity  and  then  subtracted  from  the  above  sum.  The  ignition  delay  is 
then  determined  by  multiplying  this  time  by  the  shock  wave  density  ratio  to 
convert  from  laboratory  time  to  particle  time.  The  appearance  of  the  timing 
marks  simultaneously  on  both  oscilloscopes  helps  to  considerably  reduce  data 
reduction  errors . 

The  use  of  emitted  radiation  as  a  characteristic  of  the  initiation  and 
progress  of  a  chemically  reacting  system  also  presents  several  difficulties , 
particularly  if  one  is  interested  in  radiation  from  the  OH  molecule.  At  the 
temperatures  of  interest,  almost  all  of  the  OH  is  in  the  electronic  groundstate, 
and  absorption  measurements  can  be  directly  and  quantitatively  related  to  the 
concentration  of  OH.  It  has  been  amply  demonstrated  that  below  2000°K  the 
mechanism  for  emitted  OH  radiation  is  chemiluminescent  in  nature  (Refs.  3,  4  ). 
Therefore,  before  it  is  possible  to  realte  OH  emission  to  a  specific  chemical 
composition,  it  is  necessary  to  know  the  exact  chemical  mechanism  by  which 
electronically  excited  OH  is  produced. 

The  reactions  which  are  sufficiently  energetic  to  produce  OH  in  the  ^2 
state  are  listed  in  Table  2  •  During  the  preignition  phase  of  the  hydrogen- 
oxygen  reaction,  the  concentrations  of  the  free  radicals,  O,  OH,  and  H  grow 
exponentially  with  the  form 
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References  3,5  have  shown  that  during  this  exponential  growth  period, 
the  emission  intensity  from  OH  takes  the  form  I  =  Iq  enat  where  n  is  the 
number  of  free  radicals  taking  part  in  the  excitation  reaction.  Therefore, 
comparing  values  of  the  OH  emission  time  constant  with  the  time  constant 
for  the  growth  of  any  of  the  free  radicals  will  determine  n.  The  experimental 
results  presented  in  Ref.  3,  which  compare  the  growth  constant  for  OH  emission 
with  that  of  OH  absorption,  and  Ref.  4  ,  where  the  more  sensitive  end-on 
technique  was  used  to  compare  the  growth  constant  for  OH  emission  with  the 
growth  constant  of  emission  resulting  from  the  O-CO  reaction,  indicate  that  the 
value  for  n  is  2.  Further,  the  most  probable  reaction  which  involves  2  free 
radicals  and  produces  excited  OH  was  concluded  by  both  authors  to  be 
O  +  H  +  (M)  =  OH*+  (M).  Neither  author  was  able  to  conclusively  demonstrate 
whether  this  reaction  required  a  third  body  or  not. 

The  result  of  the  more  sensitive  of  our  experiments  are  shown  in  Figure  5 
which  compares  measured  growth  constants  for  OH  emission  with  theoretical 
calculations  of  the  free  radical  growth  constant  obtained  using  the  reaction 
mechanism  given  in  Table  1.  The  factor  of  2  difference  between  our  experimental 
and  theoretical  calculations  tends  to  corroborate  the  2  free  radical  reaction 
mechanisms  proposed  in  References  3  and  4,  However,  using  this  mechanism 
we  have  been  unable  to  analytically  reproduce  qualitative  emission  intensity 
profiles  which  agree  with  those  obtained  experimentally  throughout  the  entire 
combustion  region.  This  inability  to  analytically  reproduce  the  experimental 
results  cannot  be  explained  by  introducing  a  quenching  mechanism  and 
therefore  indicate  either  a  shifting  chemical  mechanism  for  producing  excited 
OH  once  the  ignition  transients  are  passed  or  significant  third  body  effects 
which  come  into  play  as  the  reaction  proceeds. 

All  of  the  Intensity  profiles  exhibited  a  sudden  shift  in  emission  growth 
rate  at  some  time  after  Ignition.  This  shift  is  barely  perceptible  at  the  high 
temperatures  (~  2000°K)  and  becomes  quite  pronounced  as  the  temperature 
decreases.  The  appearance  of  this  "shoulder"  has  been  observed  in  OH 
emission  profiles  in  high  temperature  methane -oxygen  mixtures  under  similar 
test  conditions,  i.e.,  high  argon  dilution,  T  <  2000°K,  by  Bowman  and  Seery 
(Ref.  7).  They  proposed  that  the  chemiluminescent  reaction  producing  excited 
OH  during  hydrocarbon  oxidation  is 
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CH(sn)  +  o„(3S  )  -*  co(1F/f)  +  oh(2i:+) 

£  8 

and  the  shoulder  is  caused  by  the  competition  for  CH  between  this  reaction  and 

c(3p)  +  cH(an)  -* c,(3n  )  +  h(?  s) 

c  Q 

This  conclusion  was  supported  by  their  experiments  which  showed  a  rapid 
increase  in  C 2  starting  at  the  time  the  OH  emission  shoulder  was  reached. 

A  number  of  the  tests  in  the  present  study  were  rerun  with  the  hydrogen 
replaced  by  argon  in  an  attempt  to  determine  if  some  hydrocarbon  contaminant 
was  responsible  for  a  portion  of  the  measured  OH  emission  profile,  including 
the  shoulder.  No  emission  was  measured  in  these  tests. 

In  view  of  these  results  and  those  of  References  3  and  4  in  regard  to  the 
mechanism  of  OH  chemiluminescence  in  H2-02  systems  ,  it  appears  that  in 
the  present  case  the  "shoulder"  is  due  either  to  a  reaction  involving  O  or  H 
becoming  prominent,  or  a  complete  change  in  the  chemiluminescent  reaction. 
These  possibilities  are  currently  being  investigated. 

The  OH  emission  intensity  profiles  were  also  used  to  measure  "ignition 
delay  time"  (defined  as  the  first  appearance  of  an  emission  signal)  for  the 
hydrogen-oxygen  reaction  as  described  above.  This  data  was  used  to:  compare 
results  from  our  new  facility  with  experimental  results  obtained  by  other 
investigators;  verify  the  ignition  delay  time  predictions  of  the  computer 
analysis  with  experimentally  measured  delay  times;  and  to  provide  base  runs 
for  additive  effectiveness  experiments .  Figure  7  compares  experimentally 
measured  ignition  delay  time  for  a  H2/02  =  1  mixture  ratio  with  the  analytically 
predicted  ignition  delay  times  determined  by  the  computer  program.  The 
correlation  used  is  that  proposed  in  Ref.  6.  As  can  be  seen,  the  analytical 
predictions  are  in  excellent  agreement  with  the  experimental  measurements. 
Experiments  are  currently  being  planned  to  provide  data  on  additive  effective¬ 
ness  for  comparison  with  the  analytically  obtained  effects  of  additive  addition. 
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TABLE  1 


CHEMICAL  REACTIONS  USED  IN  COMPUTATIONS 
(REVERSE  RATE  PARAMETERS  kRev  =  AT“N  e  -B/RT 


H2 

»  H  ♦  H 

t 

A*8.00E17» 

Nal.o* 

b*o  . o  * 

02 

■  0*0 

« 

A«1.90E16t 

N*o.b  » 

R«0  »  0 

♦ 

H20 

■  OH  ♦  H 

t 

A* ] , 1 7E 17  • 

N*0.0» 

BaO.O 

* 

LEEDS 

2 

VQ 

H02 

■  02  ♦  H 

t 

A*1.59E15» 

N*0.0  » 

Ba-1.0 

» 

LEEDS 

3 

OH 

*  0  ♦  H 

• 

A*2»  OOE 18  * 

Na 1 »  U  ♦ 

BaO.O 

* 

Oh  *  OH 

■  H202 

• 

Aal,l7E17, 

NaO.O  , 

Ba«*5.5 

* 

LEEDS 

3 

Q 

ohx  *  photon  *  o* 

Ht 

Aal.OOtl^* 

N»0.0» 

BaQ.O 

t 

END  TOP 

REa* 

HgO  ♦  H 

■  Hg  ♦  OH 

• 

Aa2,l9E13» 

N*0 • 0  * 

Ba5.l5 

* 

LEEDS 

2 

X 

♦ 

X 

o 

*  H2  ♦  0 

f 

A* 1 .  74E 13* 

N*0»  0* 

8*9.45* 

LEEDS 

2 

c 

♦ 

X 

o 

*  02  ♦  H 

« 

Aa3, OOE 14  » 

N*0. Oi 

Bait, 8. 

OH  ♦  OH 

■  HgO  ♦  0 

f 

AaS.75E13» 

NaO.O  » 

B*l8.C 

* 

LEEDS 

2 

VQ 

H20  ♦  H02  •  H202 

♦  0Mf 

Aa 1 • OOE 1 3  * 

Nao.U  l 

B«1 .6 

« 

LEEDS 

3 

R 

H202  ♦ 

H  ■  H2  ♦ 

H02t 

Aa9.60E12* 

NaO.O  . 

Ha24,0 

* 

LEEDS 

3R 

OH  ♦  PHOTON  *  OHX 

t 

Aal.flEb, 

N*O.Ot 

8*0.0* 
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TABLE  2 


POSSIBLE  OH  (#T5)  PRODUCTION  REACTIONS  (*E  ~  S8  Kcal/mcle) 

H  +  H  +  OH  -«  OH*  +  H2 
O  +  O  +  OH  -»  OH*  +  02 
H  +  OH  +  OH  -*  OH*  +  H20 
O  +  H2  +  OK  ■*  OH*  +  H20 
H  +  02  +  H2  -»  OH*  +  H20 

O  +  H  (+M)  -*  OH*  (+M) 

OH*  DESTRUCTION  REACTIONS 
OH*  +  M  -*  OH  +  M 

OH*  -*  OH  +  hv 
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Figure  2.  Effect  of  Adding  H  Atoms  on  Reaction  Production  Rates 
H2/02  “  1 
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Figure  5 
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COMPUTER  PROGRAM 


APPENDIX  A 


1 .  ANALYSIS 

The  computer  program  described  in  this  report  is  based  on  analysis 
developed  previously  by  the  authors  for  the  purpose  of  predicting  delivered 
specific  impulse,  including  the  effect  of  kinetic  losses,  for  liquid  propellant 
rocket  engines.  Results  of  these  earlier  efforts  are  reported  in  references  1-3. 
This  computer  program  has  been  written  for  the  purpose  of  screening  and  analyzing 
generalized  kinetic  processes.  The  method  of  solution  consists  of  integrating 
the  conversation  equations  for  the  chemical  system  in  a  form  such  that  the 
chemistry  is  generalized  for  binary  exchange  and  dissociation-recombination 
reactions.  Solid  or  liquid  phase  products  are  not  considered. 


Conservation  Equations 

The  SCREEN  computer  program  integrates  a  set  of  simultaneous  differential 
equations  along  a  pressure  defined  streamtube  (i.e.  ,P(x)  and  dP(x)/dx  are  known). 
These  differential  equations  represent  the  conservation  of  species,  mass, 
momentum,  and  energy  for  the  system  as  expressed  by  equations  (5),  (6),  (7), 
and  (8)  below. 


The  conservation  equations  governing  the  Inviscid  flow  of  reacting 
gas  mixtures  have  been  given  by  Hirschfelder,  Curtiss,  and  Bird4,  Penner5  and 
others.  The  following  basic  assumptions  are  made  in  the  derivation  of  these 
equations . 

Mass  (m,  &),  momentum  (M),  and  energy  (H)  addition  rates 
are  defined  Tor  the  system. 

.  .  The  gas  is  Inviscid. 

.  .  Each  component  of  the  gas  is  a  perfect  gas. 

.  .  The  internal  degrees  of  freedom  of  each  component  of  the  gas 
are  in  equilibrium. 


In  one-dimensional  flow,  the  conservation  equations  have  the  form* 


species 

Mtl  +  ralcJ  ■  SI  +  <l+5)  pV 

(1) 

mass 

dk(l  +  ra)  = 

momentum 

£L<1  +  *>v]  =  m-^^ 

(3) 

energy 

db  [(1  +  ™)hT  j  "  '  hT  =  ill  ci  hi  + 

(4) 

*The  independent  variable,  x,  is  taken  as  unitless  with  r*  as  the  conversion 
tactor  to  units .  The  quantity  1  +  m  represents  the  streamtube  mass  flux  normalized 
by  the  initial  streamtube  mass  flux,  i.e.,  1  +  m  =  (oVa)/(pVa) 

o 

A-l 


if  the  expansion  process  is  specified  by  the  pressure  distribution  as  a  function 
Of  distance.  Equations  (1-4)  nan  ha  written  as 


£fi 

dx 


Sj  -  rfiCj  r* 

1  +  m  +  pV 


d V  =  M  -  ftiV  _1_  dP 
dx  1  +  m  pV  dx 


=  X 

dx  C. 


r  H  -  to 


^X-  _  V(M  -  rftVl 
1  +  m  1  +  m 


1  dP  f  ffi_  ' 
p  dx  "i*l  hi  dx  j 


d£.= 

dx 


where 


1.  dP  1  d T  1  /  A  D  d'cl  N 

P  dx  T  dx  R  V  fel  Ki  dx  )  J  p 

°p"&  cicPi'  y  =  (c^R) 

fT 

Jo 


hl  = 


CpidT  +  hlO 


For  each  component  of  the  gas,  the  equation  of  state  is 

pi =  w 

Summing  over  all  the  components  of  the  mixture,  the  over-all  equation  of 
state  is  obtained 

n 

P  =  pRT  =  pT  c1Rl 

The  net  species  production  rate  for  each  species  (component) 
is  calculated  from 

n 


where 


"i*  "V  i’=i  *VVxi 

xi"  Mi 


The  equilibrium  constant,  Kj ,  is 


Kj  =  e 


-A  FA  T 


AF  Si  fiMlj  "Si  flty'lj 


(5) 

(6) 


(7) 

(8) 

(9,10) 

(U) 

(12)  • 

(13) 


(14) 


(15) 


Where  X  depends  on  the  order  of  the  reaction  and  is  calculated  only  for 
dissociation-recombination  reactions. 


(16) 


A-2 


The  computer  program  considers  chemical  reactions  defined  by 
the  generalized  chemical  reaction  equation 

where  u..  andtA.  are  the  stoichiometric  coefficients  to  be  used  in  equation 
(15)  while  represents  the  symbol  for  the  i  chemical  species. 

The  reaction  rates,  kj,  for  the  jtk  reaction  appearing  in  equation 
(15)  are  represented  In  the  Arrhenius  form 


k,  =  a,T 


l,  (-b./RT) 
J  e  J 


where 


is  the  pre-exponential  coefficient 


n,  is  temperature  dependence  of  the  pre- 

J  exponential  factor 

bj  is  the  activation  energy 

Since  each  dissociation-recombination  reaction  has  a  distinct 
reaction  rate  associated  with  each  third  body,  the  net  production  rate  for 
each  dissociation-recombination  reaction  should  be  calculated  from 

A  T  n  yn  n  *',h  1 

XJ  “l^l  L  Kj  Slci  "  p£l  C1  J  ckkkJ  (19) 

C 

rather  than  equation  (15).  However,  Benson  and  Fueno  have  shown  theoretically 
that  the  temperature  dependence  of  recombination  rates  is  approximately  inde¬ 
pendent  of  the  third  body.  Assuming  that  the  temperature  dependency  of  recom¬ 
bination  rates  is  independent  of  the  third  body,  the  recombination  rate  associated 
tVi 

with  the  k  species  (third  body)  can  be  represented  as 

-n.  -b./RT 
kkj  =  akjT  e 

where  only  the  constants  a^  are  different  for  different  species  (third  bodies). 

From  equation  (19)  it  can  be  shown  that 


tn  JA  n  iA4“|f  n 

KJ  £l  C1  '  p  £l  cl  JU 


tlL  c  x 

ak) 

-n.  -b./RT 
akJT  Je  J 


(21) 


Thus,  the  recombination  rates  associated  with  each  third  body  can  be 
considered  as  in  equation  (15)  by  calculating  the  general  third  body  term  as 

n 

MJ  =  Si  mjicl 

where  m^  is  the  ratio  a^/a^j  • 


In  order  to  numerically  integrate  the  above  equations  (1,  5,  6,  7)  it  li 


necessary  to  input  the  following  type  of  information  concerning  the  chemical 

system: 

Boundary  Conditions: 

xo 

initial  axial  position 

Po 

initial  pressure 

To 

initial  temperature 

V 

o 

initial  velocity 

xmax 

final  axial  position 

P(x) 

table  of  pressure  vs  axial  position 

dP(x)/dx 

table  of  pressure  derivatives  vs  axial  position 

Species  Information*: 

Mi 

species  name 

mi 

species  molecular  weight 

CP1(T) 

species  specific  heat 

hj(T) 

species  enthalpy 

ft(T) 

species  free  energy 

*The  items  Cpj,  hj,  and  are  not  available  directly  in  the  appropriate  units. 
The  computer  program  calculates  these  items  for  each  species  from  the  JANAF 


data  for: 


and 


H  H  298 

P°_  13° 

f  n  298 


vs  T 
vs  T 

vs  T 


AHC 


‘298 


(22) 
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Reaction  Information: 


each  reaction  must  be  Input  In  terms  of 
Its  stoichiometric  coefficients  and 
constituent  species 

constants  defining  k.  the  reaction  rates 

)  i 

third  bod/  reaction-rate  ratios 

normalization  factor  for  x 


Mass,  Momentum,  Energy,  and  Species  Addition  Functions: 

m  (x)  mass  addition  rate,  per  unit  normalized 

length,  per  unit  Initial  streamtube  mass  flux 

M  (x)  momentum  flux,  per  unit  normalized  length 

per  unit  Initial  streamtube  mass  flux 

H  (x)  energy  addition  rate,  per  unit  normalized 

length,  per  unit  Initial  streamtube  mass  flux 

AW 

fi.  (x)  i  species  mass  addition  rate,  per  unit 

normalized  length,  per  unit  initial 
streamtube  mass  flux 


A  considerable  amount  of  data  (such  as  JANAF  tables  defining  Cpj,hj, 
and  fj  ,  reaction  rate  parameters  and  cards  defining  chemical  reactions)  are 
available  with  the  computer  program.  Details  concerning  input  to  the 
computer  program  are  given  In  Section  6. 
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NUMERICAL  METHOD 


It  has  been  shown  (e.g.  reference  3)  that  explicit  methods  of 
numerical  Integration  are  unstable  when  applied  to  relaxation  equations 
(such  as  equations  (1),  (5) ,  (6) ,  (7))  unless  the  integration  step  size  Is  of  the 
order  of  the  characteristic  relaxation  distance.  Since  In  the  near  equilibrium 
flow  regime  the  characteristic  relaxation  distance  is  typically  many  orders  of 
magnitude  smaller  than  characteristic  physical  dimensions  of  the  system  of 
interest,  the  use  of  explicit  methods  to  integrate  relaxation  equations  often 
results  in  excessively  long  computation  times.  An  implicit  integration  method 
which  Is  inherently  stable  In  all  flow  situations  (whether  near  equilibrium  or 
frozen)  is  therefore  used  by  the  computer  program.  With  this  method  step 
sizes  which  are  of  the  order  of  the  physical  dimensions  of  the  system  of 
Interest  can  be  used,  reducing  the  computation  time  per  case  several  orders 
of  magnitude  when  compared  with  conventional  explicit  integration  methods. 


Consider  N  first  order  simultaneous  differential  equations 


dyi 

-j— -=  f.  (x,  y,  ,  .  .  . ,  yM) 


1=1,  2,  . . .,  N 


dx  "  ‘i  rl  '  *  *  *'  rN' 
with  known  partial  derivatives  (l.e.  the  Jacobian  for  the  system)* 

*f. 


al  =S5T 

Ml 

The  following  implicit  difference  equations  are  used  by  the  computer  program 
to  determine  the  y^  n+i,(the  subscript  n  denotes  the  n**1  integration  step) 


l 


yl,n+l  yi,n  +  ^l,n+l 


'h=xn+rxn 


where 


‘i,n+l 


=  I"  f,  +  a.  h  +£.  k.  .  ]  »h 
[_  i,n  l,n  j=l  i,j ,n  J,n+1  J 


for  the  initial  step  and  for  restart  (1st  order), 


♦The  computer  program  uses  analytic  expressions  for  calculation  of  the 
partial  derivatives,  a., 


(27) 


ki,n+i  "  3  ^kl,n  +  2  {fi,n+  “l,nh+  $1  ^i,J,nkj  ,n+l}  ‘ 


for  equal  steps  (2nd  order  with  h  *  pr^vioi.ie  H) 

ha 


cl,n+l  "  J2h 


n+i+hn)  •  hn  _  kl#n  +  ^fl,n+0!i/nhn+l 


K 

+  Li  ft  W 
j=l  '°i,J/nKj,n  +  l 


fc^‘+h")] 


(28) 


for  unequal  steps  (2nd  order  with  h  i  previous  h) 

A  derivation  of  these  equations  is  given  in  reference  2. 

If  the  flow  is  frozen,  the  explicit  form  of  the  above  equations  can 
be  used  (c^  =  0,  jS^  =  0),  i.e.  equations  (26),  (27),  and  (28)  are  each 
reduced  from  an  NXN  system  of  linear  simultaneous  equations  to  N  explicit 
equations  (N  =  3  +  No.  of  species). 

Control  of  the  integration  step  size,  h,  is  provided  by  calculating 
estimates  for  the  truncation  error  and  comparing  these  to  an  Input  criterion,  6. 


The  step  size  Is  halved  if  for  any  1  =  1,  2 ,  . . . ,  N 

Et  >  fi 

The  step  size  is  doubled  If  for  all  1  =  1 ,  2 ,  .  . . ,  N 

E  <  * 

H  10 


where 


k,  .,-2k,  +  k.  _  , 

g  —  I  * # n  i/H— l 


ji 


3  k  -k 
5*l,n+l  *i,n 


The  above  expression  for  is  derived  in  reference  2. 


PROGRAM  II I  PUT 


2. 

This  progtarr,  was  developed  or,  the  CDC  6400  computer  using  the  FORTRAN 
IV  language.  The  program  uses  CDC  explicit  overlay  and  requires  a  minimum 
of  35K  core.  Conversion  to  another  computer  system  should  be  straight  forward 
provided  sufficient  core  storage  is  available.  Program  overlay  extends  three 
letels  deep. 

Input  to  the  computer  program  Is  divided  into  the  following  6  sections  de¬ 
scribed  below. 


$THERM0-:iamellst  input  which  controls  the  thermodynamic 
data  input 

THERMODYNAMIC  DATA  -  (optional) 

TITLE  CARD  -  also  serves  for  plot  labels 

SPECIES  CARDS  -  species  to  bo  considered  and  their  initial 
concentrations 

REACTION  CARDS  -  input  of  reactions  and  rate  data 

$PR0PEL  -  namelist  Input  for  a  specific  case 

A  card  listing  for  a  complete  Input  for  a  sample  case  Is  given  in  Table  3  at 
the  end  of  this  section. 


STHERM0 


Permits  the  generation  of  a  master  thermodynamic  file  or  the  use  of  a  tape 
file  previously  generated  master  thermodynamic  file, 


Variable  Value 

NUCHEM  =  1 


0 


MAXSP  “  Input 


LIST  =  1 

“  0 
LISTX  =  1 


-1 


=  0 

INTAPE  =  Input 


Description 

A  master  thermodynamic  file  will  be 
generated  for  this  case  on  tape  unit  4. 
'Species  and  thermodynamic  data  will 
be  read  from:  unit  INTAPE  (nomlnalE55, 
the  input  file)  and  a  new  thermodynamic 
file  will  be  generated  on  unit  4.  The 
new  thermodynamic  file  will  be  end- 
filed  and  rewound  after  generation. 

A  previously  generated  master  thermo¬ 
dynamic  file  will  be  used  for  this  case. 
The  master  thermodynamic  file  must 
be  file  TAPE4.  No  other  input  variables 
are  required  for  this  option. 

If  a  master  thermodynamic  file  is  to  be 
generated  this  variable  specifies  the 
the  number  of  species  for  the  master 
file. 

A  list  of  species  names  for  those  species 
master  file  will  be  output. 

This  output  will  be  deleted. 
Thermodynamic  functions  (CP,  H,F) 
will  be  output  for  each  species 
(one  page  per  species,  52  lines  per 
page).  Species  names ,  molecular 
weights  and  heats  of  formation  will 
also  be  printed. 

Only  a  table  of  species  names, 
molecular  weights  and  heats  of  for¬ 
mation  will  be  printed. 

The  above  output  will  bo  deleted. 

Tape  file  from  which  thermodynamic 
data  is  to  be  read  (nominal  =  5). 


$END 


THERMODYNAMIC  DATA 

For  a  NUCHEM  =  1  option  in  $THERM<2l  the  p:ogram  will  read  MAXSP 

Master  Species  cards  contatninq:  species  symbolic  identifier,  molecular 

weight,  and  AH"r  ,  and  then  read  MAXSP  sets  of  thermodynamic  data 
1  298 

(CP,  H,F)  checking  names  and  card  sequences.  The  Master  Species  cards 
must  be  sequenced  sequentially  in  columns  41-50  (110  format)  and  must 
correspond  directly  to  the  order  in  which  the  thermodynamic  data  is  to  be 
read.  Thermodynamic  functions  consist  of  10  cards  per  function,  5  values 
per  card  corresponding  to  temperature  values  of  100  •*5000°K  at  100°K  intervals. 

Table  1  lists  the  species  for  which  therm'odynamic  data  is  currently  available.  For 
species  which  do  not  appear  in  Table  1  this  input  may  be  obtained  directly 
from  the  JANAF  tables.  Table  2  is  a  sample  listing  of  the  thermodynamic  function 
input  for  the  species  Ng. 


Master  Species  Cards 

Col  Information 

1-10  Not  used 

11-16  Species  Symbolic  Identifier,  6  alphanumeric 
characters  (left  justified) 

17-20  Not  used 

21-30  Species  Molecular  weight  (F)  format 

31-40  Species  AH°_  (F)  format,  K  cals/mole 

*298 

41-50  Right  justified  sequence  number  used  for 

sequence  checking  on  input  (110  format). 

51-80  Not  used 

It  should  be  noted  that  a  species  is  identified  by  the  name  assigned 
by  the  user.  In  general  this  name  is  the  chemical  symbol  e.g.  OfC^H,^. 
However,  it  may  be  useful  to  define  a  dummy  species  with  all  the  properties 
of  another  species  but  which  may  be  treated  in  a  special  manner,  e.g. 
the  percentage  of  the  total  amount  of  a  species  which  is  designated  as  an 
inert  (possibly  to  simulate  incomplete  mixing  or  combustior).  This  may  be 
done  by  defining  species  O  and  OX  where  OX  is  identical  to  O  except  in  name, 
but  docs  not  appear  in  any  reaction. 
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Col  Information 

I  Not  used 

2-10  Function  Value  at  (100+500  (rs  —  1) )  °K,  n=  card  number 

II  Not  used 

12-20  Function  Value  at  (200+500  (n-1)) 

21  Not  used 

22-30  Function  Value  at  (300+500  (n-1))  CK 
31  Not  used 

32-40  Function  Value  at  (400+500  (n-1))  °K 
41  Not  used 

42-50  Function  Value  at  (500+500  (n-1))  °K 

51-60  Not  used 

61-66  Species  Symbolic  Identifier,  left  justified 

67-68  Not  used 

69-70  Function  Definition  CP,  H,  orF,  left  Justified 

73-76  The  word  CARD 

77-78  Card  number  1-10  right  justified 

79-80  Not  used 

Species  Symbolic  Identifier,  Function  Definition,  and  Card  Numbers 
are  checked  for  consistency  on  input. 
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TITI.r,  C.VN  ' 

The  title  card  contains  free  field  Information  which  will  be  written  as  a 
header  label  for  the  program  output.  The  first  forty  characters  will  be  written 
as  a  label  on  each  plot  if  plotting  is  requested. 

SPEC  IKS  CARDS 

This  input  is  prefixed  by  a  single  card  with  SPECIES  in  columns  1  to  7 
and  with  the  words  MASS  FRACTION'S  or  M0LE  FRACTIONS  in  columns  9-22.  If 
the  identifier  for  mass  or  mole  fractions  is  omitted,  mass  fractions  are  assumed. 

Up  to  40  species  cards  may  bo  input.  Only  those  species  specified  bv  input 
soecies  cards  will  be  considered.  The  order  of  the  Input  species  cards  is  independent 
of  the  order  in  which  the  species  appear  on  the  master  thermodynamic  data  file. 
However,  the  order  of  the  input  species  cards  does  define  the  species  order 
for  the  specific  calculation  and  other  input  referencing  individual  species  must 
refer  to  the  order  of  the  input  species  cards.  Species  Cards  are  described  below: 


Col 

Function 

1-10 

Not  used 

11-16 

Symbol  (left  justified) 

17-20 

Not  used 

21-60 

Value  of  initial  species  concentration  (if  zero 
must  be  input  as  0.0)  free  field  F  or  E  format 

61-80 

User  Identification  if  desired 

Symbols  For  Species  Identification 

A  chemical  species  is  identified  symbolically  by  6  alphanumeric 
characters  as  follows: 

The  species  symbol  must  agree  with  that  given  on  its 
master  species  card  (columns  11-16)  used  in  generating 
the  master  thermodynamic  file. 

If  the  species  is  ionized,  the  degree  of  ionization  is 
indicated  by  a  +  (or-)  sign  followed  by  an  integer 
describing  the  degree  of  ionization  (if  no  integer  is 
given  the  species  will  be  assumed  singlely  ionized.) 

The  species  symbol  may  not  contain  the  characters 
*  or  =  .  The  special  species  symbol  PII0T0N  is  reserved 
for  specifying  radiative  reactions. 
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examples: 


Symbol 

Interpretation 

CL 

Cl 

NA+ 

Na+ 

K+2 

K++ 

CL2-2 

Cl'' 

H202 

H2°2 

REACTION  CARDS 


This  input  is  prefixed  by  a  single  card  with  REACTI0NS  in  columns 
1  to  9.  Up  to  50  dissociation  reactions  and  a  total  130  reactions  may  be 
input  following  this  card.  Only  one  card  per  reaction  is  permitted.  Cards 
specifying  dissociation-recombination  reactions  must  precede  cards  spec¬ 
ifying  exchange  reactions.  The  content  and  format  of  the  reaction  cards 
is  defined  as  follows: 


Each  card  is  divided  Into  five  fields,  separated  by 
commas.  Each  field  contains: 


field  1 

the  reaction 

field  2 

A=  followed  by  the  value  of  A 

field  3 

N=  followed  by  the  value  of  N 

field  4 

B=  followed  by  the  value  of  B, 

the  activation  energy  (Kcal/mole) 

field  5 

available  for  comments 

Rules  for  specifying  the  reaction  are  given  below. 

The  values  A,  N,  and  B  define  the  reverse  reaction  rate, 
k,  as 

Jc-A.  I'N. 


All  three  reaction  rate  parameters  must  be  Input.  The 
numeric  value  of  each  parameter  may  be  specified  in 
either  I,  F,  or  E  format.  If  E  format  is  used  the  E 
must  appear  before  the  exponent. 

There  may  be  no  blanks  between  the  characters  A  and  equal  sign, 
the  N  and  equal  sign, and  the  B  and  equal  sign. 

Input  rate  parameters  are  In  units  of  cc,  °K,  mole,  sec. 


The  general  form  of  a  reaction  Is: 

N ,  *Symbolj  +N2‘;Symbol2+. . . .  =  Na*Symbola  +  Nj3*Symbolj3+.  . 

where  the  left  hand  side  represents  reactants  and  the 
right  hand  side  represents  products. 
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Each  Symbol  must  be  as  defined  on  an  input  species  card 
(see  the  description  of  SPECIES  CARDS). 

The  multipliers,  N,  must  he  Intergers  and  represent 
stoichiometric  coefficients.  If  no  stoichiometric 
coefficient  is  given,  the  value  I  is  assumed.  It 
is  required  that 

N:  +  N2  +  ....  <10 

and 

Na  +  +  . . . .  s  10 

examples: 

Reaction 

NA++CL-  =  NACL 
B+2+M-2  “  BM 
BE+2+2*0H-  -  BE0H0H 


Interpretation 
Na++  CZ  -  NaCje 

B++  +  M”  »  BM 

Be+  +  20H “  «  Be(0H)2 


A-15 


The  di:.: .Delation  -  recombination  reactions  specifying  third  body  terms 
must  precede  other  types  of  reactions,  and  must  be  followed  by  the  directive: 

Col  1  J. 

end  rmt  UKAX 

all  reactions  prior  to  the  above  directive  will  have  a  third  body  term  added  to  each 
side  of  the  inaction.  E.  g. 

H2  -  II  HI,  .  .  . 

END  TBRRCAX 

is  the  same  as 

H  +M  =  H  +  H+M,  .  .  . 

where  M  is  a  generalized  third  body.  Specific  third  body  effects  may  be  included 
by  inputing  specific  third  body  reaction  rate  ratios  XMM  (J,  I) . 

Radiative  reactions  may  be  considered  using  the  special  species 
PH0T0N.  The  PH0T0N  may  only  appear  on  the  left  hand  side  of  the  equal  sign 
for  a  reaction. 


The  reaction  set  is  terminated  by  a  card  containing 
IAST  CARD  in  columns  1  to  9. 
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$PPOPEL 


RSTAR 


NOCHEM 


Initial  normalized  axial  position 
Initial  pressure 
Initial  temperature 
Initial  velocity 

Axial  distance  normalizing  factor 
Normalized  axial  distance  for  run  termination 
-  1  if  a  frozen  chemistry  case  is  desired 


Unils. 

None 

PSIA 

°R 

Ft/Sec 

Inches 

None 

None 


ITP  =  1 

IRQ^P  -  1 

ICHEMP  (1)  = 


DCP  =  1 

MLP  -  1 

MSP  =  i 


Plotting  requested 
Plot  temperature 
Plot  density 

Species  numbers  for  desired  species 
concentration  plots,  up  to  30  species 
Plot  functions  vs  normalized  distance 
Plot  species  M0LE  fractions 
Plot  species  MASS  fractions 


Print  every  ND3rd  step  beginning  with  the  NDlst  step  until 
the  ND2nd  step  (The  initial  conditions  and  the  EXIT  point 
are  always  printed) 


Integration  Variables 
HI  Initial  normalized  step  size 

HMIN  Minimum  normalized  step  size 

HMAX  Maximum  normalized  step  size 

DEL  Relative  error  criterion 


JF  =0  all  variables  considered  for  step  size  control 

=  1  only  fluid  dynamic  variables  considered  for  step 
size  control. 


The  variable  which  controls  the  step  size,  i,  e.  has  the 
maximum  relative  error,  is  printed  in  the  normal  output  unde 
Integration  Parameters,  labeled  Governing  Equation.  The 
numberivariablc  correspondence  is  as  follows:  l=T;2-RO; 

3  =  V;  4  =  Species  1;S  =■  Species  2;  .  .  .  NSP  +  3  =  Species 
NSP. 
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NTB 

IPFLAG 


PSCALE 
ZTB  ( 1) 

PTB  ( 1) 

DPTB  (1)  » 


number  of  input  table  entries  for  pressure  table  (  101) 

Detei mines  the  type  of  differentiation  used  to  obtain 
derivatives  for  all  input  tables.  Reference  should  be 
made  to  Mass,  Momentum,  Energy  and  Species  Addition 
Functions  for  other  input  tables  controlled  by  JPFIAG. 

~  1  derivatives  of  input  table  obtained  by  simple  difference- 

formulae  with  second  derivatives  defined  as  0- 
=  2  derivatives  of  input  tables  obtained  by  SPLINE  fit 
(NTB  s  20) 

=  3  derivatives  input  along  with  tables 
=  4  derivatives  of  input  tables  obtained  by  parabolic 
differentiation  (NTB<20) 

Multiplicative  constant  for  the  input  pressure  table 
Normalized  axial  positions  for  input  tabular  values  for 
pressure  table  (always  Input) 

Pressure  table  (PSIA)  (always  Input) 

Pressure  table  derivative  (for  JPFLAG  =  3  option). 
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Miscellaneous  Variables 


XMM(J,  D  Reaction  rate  ratio  effect  on  reaction  J  of  species  I 

TU(1)  Temperature  above  which  approximate  extension  of 

JANAF  tables  will  occur  for  each  species  (nominal  «*  9000°R) 

IF1AST  For  overlay  reasons  must  be  set  =  1  on  the  last  case 

of  a  run 

TST0P  Time  at  which  a  run  will  arbitrarily  be  terminated 

(CP  time) 


.Qmnut 


IDQDX 

=  1 

Print  total  derivatives 

IDXJDX 

-  1 

Print  individual  reaction  net  production  rates 

IEQ0UT 

=  1 

Print  equilibrium  constant  and  its  temperature 
dependence  in  internal  units 

IRATE 

-  1 

Print  reaction  forward  and  reverse  rates  in 

Internal  units 

IMASDX 

■  1 

For  an  IADDF  ■  1  option  print  chemical  and 
addition  function  components  of  total  derivatives 

I0PXF 

*=  1 

Print  influence  coefficient  vectors 

I0PVAR(1) 

B 

Set  consecutive  entries  equal  to  K  for  variables  for  which 
influence  coefficient  vectors  are  to  be  output,  where 

K 

Yarjablg. 

1 

V 

2 

R  0 

3 

T 

4 

• 

e 

Species  No.  1 

• 

0 

e 

• 

NSP+3 

• 

0 

Species  No.  NSP 

For  example:  I0PVAR{1)  =  4,3,  will  output  Influence  coefficient  vectors 
for  species  number  1  and  for  temperature. 
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M.ir.s,  Mem 

Variable 

IADDF 

IXTB 

IADD0P 

XADSCL 

ADDX  (I) 
ADDMAS  (I) 
ADDE  (I) 
ADDEX  (D 

ADDM0M  (I) 

NPOINT  (D 

NSPADD 
ADDSP  (I,  J) 


entnm,  "iv'i ■  ■'*  Species  Addition  Functions 

=  1  Addition  functions  will  be  input 

=  Number  of  entries  in  addition  function  tables  40 

=  0  Addition  functions  input  via  tables 

=  1  Addition  functions  defined  via  multiplicative  factors 

See  EFACT  (I),  XMFACT  (I),  SPFACT  (I,  J)  below.  Note 
that  IXTB  factors  must  be  input. 

=  2  A  subroutine  for  addition  calculation  will  be  supplied  by 
the  user 

=  Multiplicative  scale  factor  for  all  addition  functions 

Note  that  input  of  XADSCL  as  1.0/(  .V.A)  will  provide 
automatic  normalization  (per  unit  initial  streamtube 
mass  flux)  for  input  addition  functions. 

=  Table  of  normaii  zed  axial  stations  for  all  input  addition 
functions 

=  Mass  addition  rate,  tier  unit  normalized  length,  per  unit 
initial  streamtube  mass  flux  (unitless) 

=  Energy  addition  rate,  per  unit  normalized  length,  per  unit 
initial  streamtube  mass  flux  (BTU/lb) 

=  Auxiliary  energy  addition  rate,  per  unit  normalized  length, 
per  unit  initial  streamtube  mass  flux  (BTU/lb)  which  will 
be  added  to  the  energy  addition  rate  ADDE  (I) .  Note  that 
ADDEX  (I)  is  independent  of  EFACT  (D  for  IADD0P  =  1 
option  and  will  be  added  to  ADDE  (I)  before  modification 
by  XADSCL 

=  Momentum  flux,  per  unit  normalized  length,  per  unit  initial 
streamtube  mass  flux  (ft/sec) 

=  Species  number  (for  current  case)  to  relate  species  addition 

functions  to  specific  species.  See  example  under  ADDSP  (I,  J) . 

=  Number  of  entries  in  NP0INT  (I)  table 

=  Species  mass  addition  rate  per  unit  normalized  length,  per 
unit  initial  streamtube  mass  flux  (unitless) 

1=1,.  .  .  IXTB  corresponding  to  number  of  entries 
1=1,.  .  .  NSPADD 
e.g.  ,  NSPADD  =  2 

NP0INT  =  3,5 

ADDSP  (I,  1)  corresponds  to  Species  3 
ADDSP  (I,  2)  corresponds  to  Species  5 

A-20 


EFACT  (I) 

m 

For  IADDOP  »  1,  ADDE  (I)  computed  as  ADDMAS(0)*EFAC 

XMFACT  (I) 

m 

For  IADD0P  "  1,  ADDM0M  (I)  computed  as 

ADDMAS  (I)  *XMFACT  (I) 

SPFACT  (I,T) 

t= 

For  IADD0P  =  1 ,  ADDSr  (I,  J)  computed  as 

ADDMAS  (I)  *SPFACT  (I,  ]) 

IHTCfTF 

=■•  1 

Restart  flag  indicating  that  a  case  Is  being  restarted 
and  directing  the  initial  enthalpy  to  be  HTC/TX 

n  r\ 

HT0TX 

= 

y  1 

Initial  enthalpy  if  a  case  has  been  restarted  (ft  /sec  ) 

DMASDX  (I) 

For  JPFLAG  =  3  option,derivative  of  the  mass  addition 
rate  d  (ADDMAS) /dx 

DEDXT  (I) 

■ 

For  JPFLAG  ~  3  option,  derivative  of  the  total  energy 
addition  rate  d(ADDE+ADDEX)/dx 

DM0MDX  (I) 

=s 

For  JPFLAG  *  3  option,  derivative  of  the  momentum 
addition  rate  d(ADDM0M)/dx 

DSPDX  (1,1) 

= 

For  JPFLAG  =  3  option,  derivative  of  the  species 
addition  rate  d(ADDSP)/dx 
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Generalised  Oblique  Shock  Calculation 


A  generalised  oblique  shock  calculation  is  specified  by  input  of  a  pressure 
table  containing  a  pressure  discontinuity,  and  a  pointer  designating  the  shock 
location . 


NSH0CK 


SHKBUG 

SMAXIT 

SKEPS(l) 

SKEPS(2) 


Pointer  designating  the  last  entry  in  the  pressure  table 

prior  to  the  shock 

e.g.  P2  =  PTB(NSHOCK+l) 

”  PTB(N  SHOCK) 

Entries  NStlOCK  and  NSHOCK+1  in  the  pressure  table 
must  have  the  same  axial  position. 

1.0  provides  intermediate  output 
0.0  provides  no  intermediate  output 

Maximum  number  of  iterations  during  a  generalized 
oblique  shock  calculation 

Relative  convergence  criterion  for  temperature  iteration 
during  a  generalized  oblique  shock  calculation 

Relative  convergence  criterion  for  overall  iteration 
during  a  generalized  oblique  shock  calculation 


UNITS 

None 

None 

None 

None 

None 
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Normal  Shook  Stagnation  Streamline  Calculation 


The  normal  shock  stagnation  streamline  calculation  option  performs  a  normal 
shock  calculation  from  specified  upstream  anc  downstream  velocities  and  continues 


the  calculation  as  a  velocity  defined  streamtube.  UNITS 

NSTAGV  =  1  specifies  a  normal  shock  stagnation  streamline  none 

calculation 

VEL1  “  Upstream  velocity  for  the  normal  shock  calculation  ft/sec 

VTB(l)  =  Array  defining  the  velocity  as  a  function  of  ft/sec 

normalized  axial  distance.  VTB(l)  is  defined  as 
the  downstream  velocity  for  the  normal  shock 
calculation 

NVTB  =  Number  of  entries  in  the  velocity  profile  none 

NVTB  s  10 1 
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Reaction  Scre  ening  1  >nut  VarinbUm 


If  a  reaction  screening  calculation  is  requested,  the  program  performs  a  two 
pass  calculation.  The  first  pass  utilises  the  complete  reaction  set  and  determines 
those  reactions  which  must  be  retained  to  satisfy  the  input  criteria  for  each  species 
screened.  The  second  pass  redoes  the  first  calculation  with  an  edited  reaction  set  and 


providers  a  summary 

page  conijarinij  both  calculations. 

UNITS 

ISCRF 

I,  specifies  a  reaction  screening  case  for  ISCSP  (I) 
species 

none 

’iSCSP(l) 

Species  number  for  those  species  to  be  screened  <40 

none 

EPSCR(l) 

Relative  retention  criterion  for  each  species  to  be 

screened.  Defined  as  the  maximum  change  in  mass 
fraction  relative  to  production  or  destruction  of  the  species 
per  unit  normalized  length  for  ail  reactions  involving 
the  corresponding  species.  40 

none 

ISC  BUG 

1  provides  intermediate  output  during  the  reaction 
screening  procedure. 

none 
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SPECIES  RESIDING  ON  MASTER  THERMODYNAMIC  TAPE 
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90 
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91 
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46 
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99 
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6,9  360", 

6,95/0, 

9,4613. 

6 , 99H0  , 

7,0690, 

M2 

CP 

CAPD 

1 

7,1960, 

7,3500, 

7 , 5120, 

7,  6  7;i  9, 

7 , 6  1  ‘i  0  , 

M2 

CP 

C  A r-'  D 

i 

7 . 9  4  S  o  , 

8,0610, 

8,1620. 

8,2520. 

8.3300, 

M2 

CP 

CA»D 

t 

vJ 

8  ,  39HQ  , 

8 , 4  5  8  0  » 

8  ,  512  0, 

6,5590, 

8,6010, 

M2 

CP 

CAPD 

4 

8,6300, 

8,6730* 

8,7030, 

6 , 7310, 

8,7560, 

M2 

CP 

CAPD 

5 

6,7/90, 

9  ,  8  0  0  0  , 

8,3200, 

8,8390, 

8 ,8550, 

N2 

CP 

CAPD 

fc 

8,8710*, 

9 , 6660  , 

6  ,  9  C  0  0  , 

8,9140, 

8,9270, 

M2 

CP 

CA»D 

7 

8 , 9390  , 

8,9500. 

8,9620, 

8 , 9720 , 

8,9830, 

M2 

CP 

CAPD 

{> 

6,9930, 

9 . 0  C  2  0  , 

9,0120, 

9,0210, 

9,0300, 

M2 

CP 

CAPD 

9 

9,0390, 

9,0480, 

9,3570, 

9.0660, 

9,0740, 

M2 

CP 

CAPD!  3 

-1 , 3790, 

- ; 6830  * 

,0150, 

, 7100 , 

1,4130, 

M2 

H 

CARD 

1 

2,1250, 

2 , 6 5  '53  , 

3 , 5960, 

4.3530, 

5,1290, 

M2 

H 

CARD 

2 

5,9170, 

6 , 7 1 r'  0  , 

7,5290, 

6,3500, 

9,1790, 

M2 

H 

CARD 

3 

10,0150, 

10 , 85^0, 

11,7070, 

12.5600, 

13,4180, 

M2 

H 

CAPD 

4 

14 , 2600. 

15,1460, 

16,3150, 

16,8660, 

17.7610, 

M2 

H 

CARD 

5 

18,6300', 

19,5170, 

20 , 3933, 

21,2600. 

22,1650, 

M2 

H 

CARD 

6 

23,0510, 

23.9390, 

24 , 3290, 

25,7190, 

26,6110, 

M2 

H 

CAPD 

7 

27,5050, 

28,3990  , 

29,2950, 

30.1910, 

31 ,0390, 

M2 

H 

CARD 

S 

31,9880. 

3? , 8BD0  , 

33,7880, 

34.6900, 

35,5930, 

M2 

H 

CARD 

9 

36,4960, 

37,4000, 

33, 3360, 

39,2120, 

40,1190, 

M2 

H 

C  A  R  D 1  0 

51,9570, 

46,4070. 

45,7700, 

46,0430, 

46,5610, 

M2 

r 

CARD 

1 

47,1430 , 

47,7310. 

49,3030, 

49.8530, 

49,3780, 

M2 

F 

CARD 

2 

49,6793, 

5 3 , 3570  , 

50 ,3130, 

51.2490, 

51,6653, 

M2 

F 

CARD 

3 

52,0650, 

52,4430. 

52, 9160, 

53.1710, 

53.5130, 

M2 

F 

CA»D 

4 

53,6420, 

54 ,1600. 

54,4600, 

54 , 7660, 

55,0550, 

M2 

F 

CARD 

5 

55,3350, 

55,6060. 

55,9700, 

56,1270, 

56,3760, 

M2 

F 

CARD 

6 

56,6190, 

58.35*0, 

57,0870, 

57.3120, 

57,5320, 

M2 

F 

CARD 

7 

57,7473, 

57,95/0. 

53,1620, 

58,3620, 

58.5590, 

M2 

F 

CARD 

B 

5 9 , 751  0  . 

53,9400, 

59,1240, 

59.3050, 

59,4620, 

M2 

F 

CAPD 

9 

59,6570, 

59,3270  , 

59, 9950, 

60.1600, 

6  0 , 3220  . 

M2 

F 

CARD10 

; 


[ 


i< 
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TABLE  3 


CARD  LISTING  FOR  A  SAMPLE  CASE 


H2 

•  n 

02 

•  u  l 

A«GOM 

h20 

0.  o 

0 

0*0 

H 

0*0 

03 

0*0 

OH 

0  •  0 

OHX 

0  •  0 

H02 

0. 0 

H202 

0.0 

photon 

0.0 

REACTIONS 


H2  *  H  ♦  H 

f 

A*6.00E17* 

N*  1 . 0  » 

0*0.0* 

02  ■  0  ♦  0 

f 

A*1.9l)Elbi 

NaU.b  * 

B*0  • 0 

• 

H20  ■  OH  ♦  H 

9 

Aa 1 . 1 7E 1 7 » 

N«(l. Ot 

B*  0  #  0 

9 

LEEDS 

2 

VO 

H02  «  02  ♦  H 

9 

Aal.59E15. 

NaO.O  * 

Ba-1.0 

9 

LEEDS 

3 

OH  ■  0  ♦  H 

9 

A»2 • QOE 1 8 » 

Na 1 . 0  t 

8*0.0 

9 

OH  ♦  OH  a  H202 

9 

A*1 . 1 7E 1 7  f 

Nau.D  * 

8*<*5.S 

9 

LEEDS 

3 

Q 

ohx  ♦  photon  ■  o* 

H  f 

Aa] . 00E14* 

N«0.0* 

8*0.0 

9 

END  TBR  HEAX 

H20  ♦  H  »  H2  ♦  OH 

9 

A*2  » 1 9E 1 3  » 

Na  1)  .  0  * 

8*8.15 

9 

LEEDS 

2 

OH  ♦  H  «  H2  ♦  0 

9 

Aal.74E13» 

N*  0  .  U  * 

8»9.45* 

LEEDS 

2 

OH  ♦  0  «  02  ♦  H 

9 

A-3.00E14* 

NaO  ,U  « 

8*16,8* 

OH  ♦  OH  *  H20  ♦  0 

9 

Aa5.75El3» 

N*O.U  • 

B»  18.0 

9 

LEEDS 

2 

VQ 

H20  ♦  H02  ■  H202 

♦  OH  * 

A«] ,00El3t 

N*0.y  * 

0*1.8 

9 

LEEDS 

3 

H 

H202  ♦  H  *  H2  ♦ 

HOP  » 

A*9.60E12» 

NaO.tl  » 

8*24  •  C 

t 

LEEDS 

3R 

OH  ♦  PHOTON  *  OHX 

9 

A»1 • 8E6* 

NaO . 0  * 

8*0.0* 

A-?7 


TABLE  3  (Continued) 


SPROPCl 

2*0.0*  PI". 735.  T*527.4,  V* 4385,0*  RbT Ah*  12.0*  EA11-.3* 
NU3«S* 

Hla.002*  HM I N« .  002  *  HMaX*0.2*  OEL*0.U01*  JF*0* 

NTB*4 »  JPFLA6"1*  PSCAIE*.7350*2TbU>*0.0*  0.01,  0.01*  100.0* 
PTB(1)«1.0,  1.0*  22. 1*  22,1* 

TCONO<1)«40#0«0» 

XMM*20OO*l • 

X*M(1*1>*5,U,  XMM11*4)»15.0* 

XMM  (3*11*0*25*  XMM  (3*2)*0.25*  XMM  ( 3 *  3)  *u •  U<>  * 

XMM ( A  *  1 1 *5, 0  *  xmm  14 *21 *2 • 0 •  XMM (4*4) *32.5  * 

IFLAST.i, 

TSTOP*120.0 
NShOCK*2»  SMKflUG* 1 . 0 » 

IScRF*1*  I SCSP  (  1 )  a6 »  EP5CP  H  •  *.0i»  ISCt(Ub*i« 

EXIT. 0.3, 

EXIT*. 2. 

SEND 
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APPENDIX  B 


EXPERIMENTAL  EQUIPMENT 


APPENDIX  3 


1  .  GENERAL  DESCRIPTION  OF  THE  SHOCK  TUBE  FACILITY 

The  Dynamic  Science  shock  tube  consists  of  a  driver  section  diaphragm 
holder,  driven  section,  acceleration  section  and  test  section,  supported 

on  a  9  inch  wide  steel  channel.  An  overall  view  of  the  facility  is  shown  in 
Figures  B-l  and  B-2. 

a.  Driver 

The  driver  has  an  internal  diameter  of  3.250  inches,  a  wall  thickness  of 
.500  inches,  and  a  total  length  of  5  feet.  It  is  constructed  of  type  304  stainless 
steel,  the  design  operating  pressure  being  2500  psi  with  a  safety  factor  of  5. 

For  ease  of  handling  during  diaphragm  replacement  the  driver  is  mounted 
on  roller  supports. 

b.  Diaohraam  Holder 

The  diaphragm  holder  Is  of  the  free  piston  design  and  is  shown  in  Figure 
B-3.  Any  elongation  of  the  clamping  bolts  due  to  the  pressure  load  allows  the 
piston  to  move  and  maintain  the  clamping  load  on  the  diaphragm. 

To  help  prevent  the  diaphragm  petals  from  rebounding  and  creating  addi¬ 
tional  flow  disturbances  the  internal  surface  of  the  downstream  section  of  the 
diaphragm  holder  has  been  grooved  permitting  the  petals  to  deform  Instead  of 
bouncing.  To  further  insure  clean  petalling  without  loss  of  diaphragm  material, 
the  diaphragms  are  scribed  in  the  form  of  a  90°  cross  by  a  milling  machine  with 
either  60°  or  90°  included  angle  cutters.  A  typical  diaphragm  before  and  after 
rupture  is  shown  in  Figure  B-4. 

c.  Driven  and  Acceleration  Sections . 

The  10  foot  driven  section  and  the  acceleration  section,  composed  of 
two  10  foot  sections,  are  essentially  identical  in  construction.  Both  have  wall 
thicknesses  of  .500  inches  and  inside  diameters  of  3.250  inches  with  all  in¬ 
ternal  surfaces  honed  to  a  32  micro-inch  finish  to  minimize  wail  viscous  effects. 
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Figure  B-3.  Diaphragm  Holder. 


ure  B~4 .  Shock.  Tube  Diaphragm 


Allonment  of  adjacent  sections  is  maintained  by  utilizing  spigot  and  counterbore 
joints.  A  port  is  provided  in  the  end  of  the  driven  section  nearest  to  the  main 
diaphragm  for  evacuation  and  sample  gas  loading. 

To  permit  access  to  the  second  diaphragm  station  the  driven  section  is 
supported  on  roller  mounts. 

Each  10  foot  section  contains  three  equally  spaced  instrument  ports  which 
can  accommodate  either  thin  film  heat  transfer  gauges  for  shock  wave  veloalty 
measurements,  or  pressure  transducers. 

d .  Test  Section 

The  test  section  is  an  assembly  of  two  components,  a  circular  tube  with 
a  6  inch  internal  diameter  and  a  .  375  inch  wall  thickness,  and  a  mating  rectangular 
insert.  The  rectangular  insert  (shown  in  Figure  B-5)  slides  inside  the  circular 
tube  with  its  wedge  shaped  leading  edges  protruding  into  the  end  of  the  acceleration 
section.  The  observation  windows  are  mounted  on  the  rectangular  insert  and  are 
Isolated  from  the  flow  eternal  to  the  Insert  by  the  window  clamps  shown  in  Figure 
B-6. 

The  purpose  of  this  assembly  is  to  remove  a  rectangular  gas  sample  from 
the  cylindrical  flow  in  the  acceleration  section  and  then  isolate  it  from  the  re¬ 
mainder  of  the  flow.  The  plane  flow  geometry  obtained  in  this  manner  enables  a 
simplified  reduction  of  the  experimental  data. 

A  spring  loaded  piston  relief  valve  is  fitted  into  the  end  plate  of  the  test 
section.  In  case  of  an  unexpected  overpressure  this  valve  opens  a  one  inch 
diameter  port  to  atmosphere  allowing  the  high  pressure  gas  to  be  vented  outside 
the  building. 

2  .  SHOCK  TUBE  OPERATION 

a.  Pressurization  and  Vacuum  Systems 

The  shock  tube  pressurization  and  vacuum  systems  are  shown  sche¬ 
matically  In  Figure  B-7.  A  16  Inch  Bourdon  type  pressure  gauge  is  used  to  measure 
the  driver  pressure  during  the  pressurization  process.  A  check  valve  in  the  line 
connecting  the  gauge  to  the  driver  section  maintains  the  gauge  at  the  pressure 
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Figure  B-7.  Vacuum  and  Pressurization  System  Schematic. 


reached  at  the  time  of  diaphragm  rupture.  Tnls  permits  a  more  accurate  calibration 
of  dilvei  pleasure  versus  diaphragm  thickness  and  scribe  depth. 

Two  vacuum  systems  composed  of  roughing  pumps  and  2  inch  oil  vapor 
diffusion  pumps  are  used  to  evacuate  the  tube  prior  to  the  loading  of  the  sample 
gas.  One  set  of  pumps  evacuates  the  acceleration  and  test  sections  and  the 
other  evacuates  the  driver  and  driven  sections.  When  evacuated,  the  tube 
pressure  and  leak  rate  are  measured  by  Edwards  High  Vacuum  Model  M6A  Pirani 
Gauge  heads. 

The  controls  operating  the  remote  actuated  valves,  and  the  Pirani  gauge 
head  read  out,  are  located  in  the  console  shown  in  Figure  B-8. 

The  sample  gas  mixture  is  prepared  by  admitting  the  correct  partial  pressure 
of  the  individual  components  of  the  desired  mixture  into  the  mixing  chamber.  The 
partial  pressures  and  the  final  total  pressure  are  measured  with  a  mercury  manom¬ 
eter. 


3 .  INSTRUMENTATION 

a .  Shock  Speed  Measurement 

The  shock  wave  velocity  is  determined  from  the  time  required  by  the  shock 
wave  to  pass  thin  film  heat  transfer  gauges  locatod  at  known  intervals  along  the 
tube.  Signals  from  these  gauges  are  suitably  amplified,  pseudo-differentlated, 
and  then  displayed  on  a  rastered  sweep  by  a  Tektronix  type  556  oscilloscope. 

The  rastered  sweep  is  obtained  from  the  method  described  in  Reference  B-l.  A 
Tektronix  type  184  Time  Mark  Generator  is  used  to  intensity  modulate  the  sweep, 
providing  an  accurate  time  base  for  the  measurement  of  the  inter-signal  time 
intervals. 

The  heat  transfer  gauges  were  constructed  by  bonding  a  thin  film  of 
platinum  on  5  mm  diameter  pyrex  rods,  using  Liquid  Bright  Platinum  No.  05 
(Hanovia  Liquid  Gold  Division,  Engelhard  Industries),  and  following  the  pro¬ 
cedures  given  In  Reference  B-2. 
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Figure  B-8.  Operating  Console 


b .  Spectroscopic  Instruments 


Thr>  ^  n  n  l  t/ti  p/1  I  incttrnmmitt’  pi  irmn  hi  ir  4  m,  ii"r\  nfifU  tV\  p  nU«/>l»  f  hUa  U-it^ 

the  capability  of  determining  chemical  specie  concentrations  by  absorption 
or  emission  spectroscopy.  The  instrumentation  (shown  in  Figure  B-9)  consists 
of  a  Warner  and  Swasoy  Model  21  Focused  Radiation  Source  and  Model  30 
Optical  Transfer  System  which  illuminate  the  entrance  slit  of  a  Hilger-Engis 
Model  600  Monochromator-Spectrometer  with  radiation  emanating  from  the  test 
section  of  the  shock  tube.  The  Monochromator-Spectrometer  is  constructed 
so  as  to  allow  the  intensity  of  several  spectral  lines  with  significantly 
different  wave  lengths  to  be  monitored  simultaneously  by  photomultiplier  tubes. 
The  photomultiplier  used  ir.  the  experiments  discussed  above  is  an  EMI  9601B 
(S-11C  photocathode)  with  the  operating  voltage  supplied  by  a  Sorenson  power 
supply . 

c .  Monochromator  Dispersion  and  Operating  Bandwidth 

The  general  equation  which  relates  the  dispersion  angle  at  a  given  wave¬ 
length  to  the  configuration  of  the  instrument  and  the  properties  on  the  grating 
used  is 

s\  =  ~  (sin  i  +  sin  cp)  (1) 

s  =  the  order  used 

3  =  the  number  of  rulings  per  °A  on  the  grating 
i  =  angle  of  incidence 
v  =  angle  of  diffraction 
For  this  instrument 

's  -  w  -  § 
i  =  w  +  $ 

Therefore,  (1)  becomes 

sX  =  \  sin  w  cos  i  w  =  sin  *  .  ?  ^  ^  .  (2) 

?  v  2  cos  $.  '  1 
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The  dispersion  in  °A  per  degree  is  found  by  differentiating  (1) 


dX  _  cos  w 
dcp  3  s 

Linear  distance  at  the  exit  focal  plane  is  related  to  the  diffraction  angle  by 

■£  ~  f  tan  cp  (3) 

f  =  focal  length  of  instrument 

Differentiating  (3)  gives 

den  cos  cn 

d  4  f 

Therefore,  the  dispersion  at  the  exit  focal  plane  in  °A/mm  is  given  by 

dX  _  _d2.  =  dcs 
di '  dc  d£ 

,,  2 

di  3  sf 

when  f  is  in  units  of  mm. 

For  the  present  instrument 

^  =  27.7  °A/mm 
at  X  «  3080  0  A  . 

Figure  B-10  presents  the  measured  intensity  of  the  mercury  2967.3  °A 
line  versus  monochromator  wavelength  setting  at  1  °A  intervals.  Using  the 
calculated  dispersion,  the  bandwidth  of  uniform  intensity  transmission  is  found 
to  be  21 .8  °A  which  agrees  very  well  with  that  obtained  from  Figure  B-10. 
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intons  it  y 


Wavelength 
£  X  =  1  CA 


Figure  8-10.  Transmission  Bandwidth  for  Mercury 
2967  °A  Line 
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Analytical  and  experimental  investigations  of  the  ignition  and  combustion  of  hydroger 
oxygen-argon  mixtures  are  presented.  A  one-dimensional  kinetics  program  with 
generalized  chemistry  and  provisions  for  mass  addition,  momentum  addition,  heat 
loss,  mixing,  shock  waves,  and  a  rate  screening  option  has  been  developed  and 
used  to  analyze  the  effect  of  free  radical  additives  on  the  ignition  delay  time  in 
hydrogen-oxygen  mixtures.  This  computer  program  has  also  been  used  to  reduce 
shock  tube  measurements  of  hydrogen-oxygen  ignition  delay  and  OH2S-3II 
emission.  The  experimental  results  indicate  that  the  reaction  O+H  (+  M)-*OH*(+  M) 
is  responsible  for  producing  OH (3£)  during  the  induction  period;  however,  this 
mechanism,  when  input  into  the  computer  program,  was  not  sufficient  to 
qualitatively  reproduce  OH  emission  intensity  profiles  obtained  experimentally. 

Work  aimed  at  determining  the  mechanism  for  OH  emission  throughout  the  reaction 
zone  is  continuing. 
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